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ABSTRACT 

We present a large-scale view of the magnetic field in the central 2° x 2° region 
of our Galaxy The polarization of point sources has been measured in the J, 

H, and Kg bands using the near-infrared polarimetric camera SIRPOL on the 

I. 4 m telescope IRSF. Comparing the Stokes parameters between high extinction 
stars and relatively low extinction ones, we obtain polarization originating from 
magnetically aligned dust grains in the central few-hundred pc of our Galaxy. 
We find that near the Galactic plane, the magnetic field is almost parallel to the 
Galactic plane (i.e., toroidal configuration) but at high Galactic latitudes (| b \> 
0?4), the field is nearly perpendicular to the plane (i.e., poloidal configuration). 
This is the first detection of a smooth transition of the large-scale magnetic field 
configuration in this region. 
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INTRODUCTION 



It has long been a subject of controversy whether the magnetic field (MF) in the central 
region of our Galaxy has a dipolar or toroidal geometry on a large scale. Previous observa- 
tions have been mainly based on sub-millimeter polarimetric observations of dense molecular 
clouds, or radio observations of non-thermal radio filaments. The former show that the MFs 
in th e dense molecular clouds run almost parallel to the Galactic plane, indicating a toroidal 
MF jNovak et al.ll2003l : Ichuss et alil2003h. By contrast, the radio filaments are aligned n early 



perpendicular to the Galactic plane (lYusef-Zadeh et al 



are h ighly polarized along the filaments' long axes (ITsuboi et al 



1984 



2004: 



LaRosa et al 



1986 



2004), and 



Yusef-Zadeh fc Morris 



19871 ). suggesting a large-scale, pervasive poloidal MF in the intercloud medium. 



The MFs both in dense molecular clouds and in the diffuse interstellar (intercloud) 
medium are observable via optical and near-infrared (NIR) polarimetry of intrinsically un- 
polarized stars with the polarization attributed to differential extinction by dust grain s 



aligned with their short axis along the direction of the MF (e.g., iMathewson fc Fordlll97CH ) . 
The polarization vectors of such stars trace the orientation of the MF projected onto the 
plane of the sky, and measurements of stars at different distances reveal the change of the 
projected MF orientation along the line of sight. Although there are numerous stars near 
the Galactic center (GC), high extinction by the intervening dust grains means that infrared 



polarimetry is ideal (e.g., iTamura et al.lll987l ) in understanding the MF configuration in the 



number of stars so far (e.g., 


Houeh et al. 


Eckart et al. 


1995; 


Ott et al. 


1999]). 



Hough et all Il978l : Kobavashi et all Il983l : Bailev et al.l Il984 



In this letter, we report obs ervations of the GC reg ion with a wide-field NIR polarimeter. 
In the first paper of this series (INishiyama et al.ll2009i ). we presented the MF configuration 
in the 20' x 20' region of the GC. In the following observations from 2007, more than 300,000 
stars are detected in the central 2° x 2° region with an error in the degree of polarization 
less than 1%. The statistical treatment of the large number of stars enables us to obtain the 
first large-scale view of the MF in the GC. 



2. OBSERVATIONS AND DATA REDUCTION 



Polarimetry for a large number of stars in the central 2° x 2° region of the Galaxy 
was obtained with IRSF/SIRPOL between July 2006 and August 2009. SIRPOL consists 
of a single-beam po larimeter (a half-wave plate rotator unit and a fixed wire-grid polarizer; 
Kandori et al.ll2006l ) and the NIR imaging camera SIRIUS (Simultaneous Infrared Imager for 
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Unbiased Survey; iNagashima et al.lll999t iNagayama et al.ll2003l ). and is attached to the 1.4- 



m telescope IRSF (Infrared Survey Facility). The detectors are three 1024 x 1024 HgCdTe 
arrays, with a scale of 0'.'45 pixel -1 . SIRPOL provides images of a 7f7 x 7'7 area of sky in 
three NIR bands, J (1.25//m), H (1.63/im), and K$ (2.14/xm), simultane ously. The filte r 



system of IRSF/SIRPOL is similar to the Mauna Kea Observatory system (jTokunagal 120001 ) . 



10-s exposures were made, at each 4 wave plate angles (0?0, 22?5, 45?0, 67?5), and at 
each of 10 dithered positions. The weather condition was photometric, with a mean seeing 
of ~1'.'3 (J), ~1'.'2 (H), and ~1'/1 (Kg). Twilight flats were obtained before and after the 
observations. The IRAF (Image Reduction and Analysis FacilityjQ software package was 
used to perform dark- and flat-field corrections, followed by sky background estimation and 
subtraction. 



ANALYSIS 



To obtain the photometric magnitudes and errors, we used the DAOPHOT package 



(jStetsoru 119871 ) . The sources detected in Stokes I images [I 



( J 0° + J 22?5 + J 45 C 



+ W/2] 



were input to the ALLSTAR task for PSF-fitting photometry. About 10 sources were used 
to construct the PSF in each image. Photometric calibrati on was performed wi th reference 
to the magnitudes of point sources in the 2MASS catalog (ISkrutskie et al.l 120061 ) . 



The Stokes parameters /, Q, and U for point sources were determined from aperture 
polarimetry, which gives a better photometric result than PSF-fitting photometry in our GC 
data. We obtained the intensity for each wave plate angle (I Q °, J 22 ° 5 , J 45 °, I 6 f 5 ) with IRAF 
DAOFIND and APPHOT tasks. The size of the aperture is slightly different among fields 
due to variations of the seeing, and an aperture size of 1.5 x FWHM was adopted, which was 
measured on Stokes I images. Our data analysis results in ~ 140, 000( J), ~ 360, 000(H), and 
~ 380,000(1^5) sources in the observed area with an error 5P in the degree of polarization 
P less than 1%, where 

p = v/(Q2 + fpj/I, 



and 



5P 




Q 



+ 



2 r 



U 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Polarimetric accuracy was checked with duplicate sources in overlapping regions of adjacent 
fields observed under different atmospheric condition. The means of the average differences 
of Q/I and U/I were 0.26 % and 0.27 % in the H band, and 0.25 % and 0.26 % in the K s 
band, respectively. 

To derive the MF configuration at the GC, we separated stars into two subgroups ac- 
cording to their colors: "blue" and "red" stars. Most of the stars detected in the observations 
are late-type giants neat the GC, and their intrinsic H - K s colors are almost the samfi 
So their observed colors depend on the amount of the interstellar extinction which in turn 
depends generally on the distance to the stars. The "red" ( "blue" ) stars suffer strong (weak) 
extinction, and thus they are distributed at the far (near) side of the GC. The light from 
stars at the near side of the GC (i.e., "blue" stars) is transmitted through the interstellar 
dust in the Galactic disk, and the polarization is produced there, whereas the light from stars 
at the far side of the GC (i.e., "red" stars) is transmitted through the dust in the disk and 
the GC. By comparing the Stokes parameters (Q/ I and U/I) of the "blu e" and "red" stars, 



we obtain polarization of the GC component (see iNishiyama et al.ll2009l . for more detail). 



As a first step in deriving the MF configuration at the GC, the field was divided into 
sub-fields of 2' x 2', and H — Ks histograms were made for each sub- field. These histograms 
were fitted with the Gaussian function, and means {mjj-Ks) an d sigmas (<th-k s ) were cal- 
culated for each sub-field. Using the mean and sigma values, the stars were divided into two 
subgroups: "blue" stars with the color of (m H _ Ks — 2cr H _ Ks ) < (H — K s ) < m H _ Ks , and 
"red" stars with m H _ Ks < (H — Kg) < {tuh-Ks + ^ a H-K s )- 

Mean {mn-Ks) an d sigma (cfh-k s ) in the H — Ks histograms were determined for stars 
with H < 15.5, which are brig ht enough to avoid the influence of the limiting magnitude, 
H ~ 16.6 (see also Fig. 7 of Nishiyama et alibood ). This means that our observations 



are equally sensitive to stars on the far and near sides of the GC, and using this criterion, 
stars on both sides of the GC can be distinguished. The "blue" and "red" stars are in fact 
distributed evenly in the 2' x 2' sub-fields. Hence, it can be safely assumed that the peaks 
in the H — Ks histograms correspond to a real peak of spatial distribution of stars along the 
line of sight. 

Next, Q/I and U/I histograms in the Ks band were constructed for the "blue" and 
"red" stars with 5P < 3% in each sub-field. Their means, < Q/I >b, < Q/I >r, < U/I >b, 



2 The standard deviation of the in trinsic H — K$ colors o f the observed giants in the GC is calculated 
to be 0.06 from the Galactic model by lWainscoat et al.l (|l992l ). This is small enough compared to the mean 
of the color ranges of red and blue stars, 2(jh-Ks — 0.38. Therefore we proceed on the supposition that the 
observed colors of the giants depend on the distance to them. 
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and < U/I >r, were used to obtain the degree P and position angle 9 of polarization for 



"red" minus "blue" (i.e., GC component) with the following equations (lGoodrichlll986[ ): 



Pi 




R-B 



£r-b = ~ arctan " ® 

The errors of < Q/I > and < U/I > were calculated from the standard error on the mean 
a/v^iV of the Q/I and U/I histograms, where a is the standard deviation and N is the 
number of stars. A few sub-fields with strong extinction are excluded from our analysis. We 
confirmed that the ratio of the number of the "red" and "blue" stars in each sub-field has a 
mean of ~ 1 even near the Galactic plane (I < 0?2), suggesting no bias effect to detect more 
stars at the near side than the far side of the GC due to strong extinction. Finally, < Q/I > 
and < U/I > are averaged in a circle of 2'A radius with a 3'.0 grid, and -Pr_b and #r_b are 
calculated with the equations (JTJ and (|2J). The resultant mean and rms of Pr_b are 0.29 % 
and 0.12 % in the Ks band, respectively. 

The strong concentration of stars towards the GC allows measurement of the MF struc- 
ture in the GC region. Most of the stars detected in our observations have an intrinsic color 



of re d giants with spectral type of K and M. We used a Galactic model (IWainscoat et al. 



19921 ) to show the degree of central concentration of the stars along the line of sight, and 
found that more than 85 % of stars within the central 2° x 2° region are expected to be within 
300 pc of the GC. So the observed position angles of polarization of the GC component show 
those of MFs, averaged along the lines of sight, in the central few-hundred pc. 



4. RESULTS AND DISCUSSION 



Fig. [TJ shows the MF direction (E- vectors of polarization of the GC component in the 
Ks band) in the GC. This is the first large-scale view of the MF in this region. The MF 
directions in the dense molecular clouds, measured with sub-millimeter (sub-mm) and far- 
infrared (FIR) polarimetry, are overplotted on the Ks data in Fig. |2j The panel (a) show s 
the MF directions derived from polarimetry in a large survey at 450 /im (INovak et al.ll2003l ). 
Near the Sagittarius A region, sub-mm polarimetry has been carried out wit h better spatial 
resolution (Fig. l21b: lDotson et al.lboool : INovak et al.lboool : Ichuss et al.ll2003h . In both cases, 
the MF directions derived from the sub-mm/FIR polarimetry are well matched with the NIR 



-6 - 



observations. It follows that the position angle of polarization derived from NIR polarimetry 
can be used to determine the direction of the MF in the GC. 

The MF direction has strong dependence on the Galactic latitude, i.e., distance from 
the Galactic plane, but not on the Galactic longitude (Fig. [3]). The MF direction (i.e., the 
position angle of polarization) is defined as the angular offset in degrees from north through 
east in Galactic coordinates. The histogram of the MF direction at | b \ < 0°A has a clear peak 
of ~ 90°, suggesting a toroidal MF configuration in this region. At higher Galactic latitude 
(| b |> 0?4), the mean MF direction appears to swing from 90° to ~ 170° which is nearly 
perpendicular to the Galactic plane. These results suggest a transition of the large-scale MF 
configuration from toroidal to poloidal in the GC region, at | b |~ 0?4. 

In the central few hundred pc region, ther e are some structures surrounding the GC, such 



as the 180-pc (or expanding) molecular ring (IKaifu et al. 



1972; 



Scovilld 119721 ). and the GC 



molecular arms composing a ring of radius of ~ 120 pc ( Sofuel 19951 ). It is therefore difficult 



to completely exclude the possibility that we are preferentially probing magnetic fields in 
these molecular structures. However, it seems to be convincing that we have detected the 
magnetic fields in molecular clouds m uch closer to the GC, from the following points. As 
shown in Fig. [2] (see also Fig. 10 of iNishiyama et al.ll2009l ). we are able to replicate the 
detailed structure of the sub-mm/FIR magnetic field maps in the circumnuclear disk, and 
two molecular clouds M-0.02 - 0.07 (+50 km s" 1 cloud) and M-0.13 - 0.08 (+20 km s" 1 
cloud). These two molecular clouds are located within several tens of p c from the GC, an d 
the possibility of contribution of foreground emission was also ruled out ( INovak et al.ll2000l ). 
Hence we can probe the magnetic fields very close to the GC. 

Previously, the MF configuration of the GC has been viewed as poloidal in the diffuse, 
interstellar (intercloud) medi um, and appro ximately parallel to the Galactic plane only in the 
dense molecular clouds (e.g., iFerrierd 120091 ) . The new data presented here shows a toroidal 
MF prevails at | b |< 0?4, eve n outside the dense mo lecular clouds (Fig. [2]), which radio 
(jOka et al.lll998f ) and sub-mm (IPierce-Price et al.ll2000l ) surveys show are mostly confined to 
within ~ 0?2 of the Galactic plane. Toward higher Galactic latitudes (| b |> 0?4), the field 
changes from a toroidal to a poloidal configuration. 

Since SIRPOL also provides images in the H band, it is possible to compare these 
with the Ks-b&nd results (Fig. [TJ). Most of the position angles of polarization in the H 
band show the same direction as that in the Kg bands within their errors. In addition, 
there is a good correlation between the degree of polarization P in the H and Kg bands 
with Pk s = 0.611 x P H . This corresponds to /3 = 1.85 for a power-law approximation 
Pa oc A _/3 which represents the wavelength dependence of interstellar polariz ation at NIR 
wavelengths, consistent with previous results toward the GC (/3 ~ 1.8 — 2.0; iNagata et al. 
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1994 ; lHatano et al.ll2010l ). The correlation coefficients in the degree and position angle of 
polarization are 0.71 and 0.92, respectively, showing strong correlation between the H and 
Ks bands. These results strongly suggest that the differential extinction by magnetically 
aligned dust grains is responsible for the observed NIR polarization. 

The existence of the toroidal MF even in the intercloud medium near the Galactic 
plane is consistent with t he observations of Galactic center diffuse X-ray (GCDX) emis- 
sion ( IKoyama et al.l I1989I ). The spatial distribution of the 6.7 keV line emission from 
Heliu m-like ions of iron is e lliptical with a major axis approximately along the Galactic 
plane (jYamauchi et al.lll990l ). Although the origin of the GCDX is as yet unclear, recent 
observations with the Suzaku satellite found a smooth and monotonic variation of the ioniza- 
tion temperature along the Galactic longitude, and a clear difference of distribution between 
the 6.7 keV line f lux and an integrate d point source flux, suggesting hot, diffuse interstellar 
plasma in origin (IKoyama et al.ll2007l ). If a poloidal MF is pervasive in this region, since the 
gravitational well of the GC region is far too shallow to confine them, the hot plasma must 
be rushing out of the plane vertically as a galactic wind. Such a vertical extension of the 
GCDX has not been observed. 

The large-scale MF configuration has important implication s for understanding the ori- 
gin of a variety of magnetic phenomena. The double helix nebula f lMorris et al.ll2006l ) extends 
almost perpendicular to the Galactic plane from b = 0?5 to 0?8. It exists in the region of 
the large-scale poloidal MF, so the nebula probably results from a dynamically ordered in- 
terstellar phenomenon involving its local interstellar MFs. On the other hand, a large part 
of the non-thermal radio filaments exist at I b |< 0?4 where the toroidal field is dominant 



(lYusef-Zadeh et al.ll2004l ; lLaRosa et al.ll2004j ). Hence the filaments do not run along the MF 
directions we obtained; rather, they are almost perpendicular to the local MF. This implies 
that it is important to reconsider the simple picture in which the non -thermal filaments are 
illuminated flux tubes within a uniform poloidal MF ( iMorrij 119901 ). for the origin of the 
filaments. 
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Fig. 1. — NIR mosaic image of the Galactic center region covering 2?0 x 2?0 in the Galactic 
coordinate, taken with the IRSF telescope and NIR camera SIRIUS. The three NIR bands 
are J (blue, 1.25/im), H (green, 1.63/xm), and K s (red, 2.14/xm). The central-parsec star 
cluster of the Galactic center is the bright yellow blob in the center of the image. Observed 
E- vectors of polarization for the Galactic center components in the Kg band are also plotted 
with cyan bars whose length indicates the degree of polarization. The vectors are averaged 
in a circle of 2'A radius with a 3f0 grid, and plotted with thick bars (detected with more 
than 3<t) and thin bars (detected with 2 — 3cr). 




Fig. 2.— Maps of the 12 C0 J— 1-0 line emission integrated over the velocity range 
Vlsr — —100 to +100 km s -1 in Galactic coordinate (jOka et al.l 119981 ). Pink bars in the 
panel (a ) are drawn parall el to the magnetic field direction derived from polarimetry at 

20031 ). and those in the panel (b) f rom measurements at 60 /im, 100 



450 urn (|Novak et al. 



/im (IDotson et al.ll2000[ ). and 350 /mi ( iNovak et al.l 120001 ) . The MF directions measured at 
near-infrared wavelength (same as Fig. [TJ) are also plotted with cyan bars whose length in- 
dicates the degree of polarization. Note that since emissive polarization in the sub-mm/FIR 
wavelengths is perpendicular to the magnetic field, vectors orthogonal to it are shown here. 
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Fig. 3. — Directions of the magnetic fields in the Galactic center as a function of the Galactic 
longitude (a) and latitude (b). The number of field per bin is shown in different colors, as 
indicated in the color scale. The MF direction is defined as the angular offset in degrees from 
north through east in Galactic coordinates. The magnetic field swings from the direction 
along the Galactic plane at b ~ 0° to the direction nearly perpendicular to the plane at 
| b |~ 1?0. Histograms of the magnetic field direction at | b |< 0?4 (red), at | b |> 0?4 (blue), 
and both (black) are also shown in the panel (c). The red histogram has a clear peak at 
the direction almost parallel to the Galactic plane, while the blue one has a peak at ~ 170°, 
almost perpendicular to the plane. 



